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Abstract 

Background: Bone morphogenetic proteins (BMPs) are members of the TGF-beta superfamily of growth factors. 
They are known for their roles in regulation of osteogenesis and developmental processes and, in recent years, 
evidence has accumulated of their crucial functions in tumor biology. BMP4 and BMP7, in particular, have been 
implicated in breast cancer. However, little is known about BMP target genes in the context of tumor. We explored 
the effects of BMP4 and BMP7 treatment on global gene transcription in seven breast cancer cell lines during a 6- 
point time series, using a whole-genome oligo microarray. Data analysis included hierarchical clustering of 
differentially expressed genes, gene ontology enrichment analyses and model based clustering of temporal data. 

Results: Both ligands had a strong effect on gene expression, although the response to BMP4 treatment was more 
pronounced. The cellular functions most strongly affected by BMP signaling were regulation of transcription and 
development. The observed transcriptional response, as well as its functional outcome, followed a temporal 
sequence, with regulation of gene expression and signal transduction leading to changes in metabolism and cell 
proliferation. Hierarchical clustering revealed distinct differences in the response of individual cell lines to BMPs, but 
also highlighted a synexpression group of genes for both ligands. Interestingly, the majority of the genes within 
these synexpression groups were shared by the two ligands, probably representing the core molecular responses 
common to BMP4 and BMP7 signaling pathways. 

Conclusions: All in all, we show that BMP signaling has a remarkable effect on gene transcription in breast cancer 
cells and that the functions affected follow a logical temporal pattern. Our results also uncover components of the 
common cellular transcriptional response to BMP4 and BMP7. Most importantly, this study provides a list of 
potential novel BMP target genes relevant in breast cancer. 
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Background 

Bone morphogenetic proteins (BMPs) are extracellular 
ligand molecules that belong to the transforming growth 
factor p (TGF-P) superfamily. To date, 21 members of 
the human BMP family have been identified [1]. BMPs 
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regulate transcription of target genes by signaling 
through type I and II transmembrane serine-threonine 
receptors. Binding of the ligand to the type II receptor 
elicits phosphorylation of the type I receptor, which, as 
a result, is able to phosphorylate other molecules and 
transmit the signal. In the canonical BMP pathway, the 
type I receptor phosphorylates receptor-regulated 
SMAD (homologue of Drosophila Mothers Against Dec- 
apentaplegic) proteins (R-SMADs, SMAD-1/5/8), which 
then bind to the common mediator SMAD4; the 
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resulting SMAD complex translocates to the nucleus to 
regulate transcription of target genes [1]. The signals 
generated by BMPs in the cell membrane may be also 
transferred into the cell via ERK, JNK and p38 mitogen- 
activated protein kinases (MAPK) [2,3]. Moreover, there 
is crosstalk between BMP signaling and other cellular 
signaling cascades, such as the Wnt, JAK/STAT and 
Notch pathways [4-6]. 

BMPs were first identified as inducers of ectopic bone 
formation in vivo [7] but were later found to be crucial 
multifunctional regulators of development [8]. During 
the last decade, the role of BMPs in cancer development 
has gained increasing interest [9-11]. The importance of 
BMP4 and BMP7 in breast cancer was highlighted in a 
survey of seven BMPs: these two ligands had the highest 
expression levels and were the most frequently 
expressed among 22 cell lines and 39 primary tumor 
samples [12]. The expression of BMP4 and BMP7 in 
breast cancer also has been demonstrated in several 
other reports [13-17]. Interestingly, BMP7 protein 
expression in primary breast tumors has been associated 
with accelerated bone metastasis formation and served 
as an independent prognostic factor for early bone 
metastasis in a study based on a set of 409 patient sam- 
ples [15] though, with a smaller set of 67 patient sam- 
ples, this association was not established [18]. 

The functional significance of BMP4 and BMP7 in 
breast cancer has been studied predominantly through 
the use of in vitro models. BMP4 was shown to inhibit 
cell proliferation in a panel of breast cancer cell lines by 
inducing a Gl cell cycle arrest [14]. The effects of exo- 
genous BMP4 on breast cancer cell migration and inva- 
sion have also been studied. For the most part, the data 
suggest promotion of these cellular abilities by BMP4 in 
several breast cancer cell lines and in normal breast 
epithelial cells [14,19], while a study in which only 
MDA-MB-231 cells were analyzed reported the opposite 
phenotype [20]. For BMP7, the results from different 
reports and different cell lines are more variable. In 
vitro examination of BMP7 manipulation have revealed 
cell line-specific effects on cell proliferation, migration 
and invasion; BMP7 induces all of these parameters in 
MDA-MB-231 cells and inhibits cellular proliferation in 
several other cell lines [21]. In opposition, in an in vivo 
xenograft mouse model of MDA-MB-231 cells, BMP7 
reduced tumor growth as well as the formation and 
growth of bone metastases [18]. 

In spite of the many years since the discovery of BMPs 
and being currently a very active topic in cancer 
research, little is known about their target genes in 
tumor conditions. The present study was designed to 
gain knowledge in this topic, by exploring the effects of 
BMP4 and BMP7 signaling on gene transcription in 
seven breast cancer cell lines and throughout a 6-point 



time series, using a genome-wide approach. We charac- 
terized the transcriptional response of breast cancer 
cells to BMP signaling in an analysis that included a 
temporal dimension and the comparison of different cell 
lines and two BMP ligands. Finally and most impor- 
tantly, we report novel potential BMP target genes rele- 
vant in breast cancer. 

Methods 

Breast cancer cell lines 

Seven breast cancer cell lines (HCC1954, MDA-MB-361, 
ZR-75-30, HCC1419, SK-BR-3, MDA-MB-231 and T- 
47D) were purchased from the American Type Culture 
Collection (ATCC, Manassas, VA) and cultured accord- 
ing to the recommended conditions except for MDA- 
MB-231 and T-47D, for which the concentration of FBS 
in culture media were 1% and 5%, respectively. 

BMP4 and BMP7 treatments 

Recombinant human BMP4 and BMP7 proteins were 
purchased from R&D Systems (Minneapolis, MN). 
Three cell lines (HCC1954, MDA-MB-361 and ZR-75- 
30) were treated with both BMP4 (100 ng/ml) and 
BMP7 (50 ng/ml) separately. HCC1419 and SK-BR-3 
cell lines received only BMP4 treatment (100 ng/ml), 
while MDA-MB-231 and T-47D were treated only with 
BMP7 (50 ng/ml). Cells were seeded on 24-well plates, 
allowed to adhere for 24 h, and treated with the BMP 
ligand or vehicle for 30 min, 1 h, 3 h, 6 h, 12 h and 24 
h (Figure 1A). Experiments were performed in triplicate 
and collected cells were pooled. 

Microarrays 

Total RNA was extracted using the RNeasy Mini Kit 
(Qiagen, Valencia, CA) and the quality of RNA was vali- 
dated using the Agilent RNA 6000 Nano Kit (Agilent 
Technologies, Palo Alto, CA, USA). Total RNA (500 ng) 
was used to generate fluorescent Cy-3-(vehicle treated 
cells) or Cy-5-labeled cRNA (BMP4- or BMP7-treated 
cells) using the Agilent Low RNA Input Fluorescence 
Linear Amplification Kit (Agilent Technologies). The 
labeled cRNAs were hybridized to the 44 K Whole 
Human Genome oligo microarrays (Agilent Technolo- 
gies) according to the manufacturer's protocol. Microar- 
ray slides were scanned (Agilent Microarray Scanner) 
after hybridization, and data were extracted using the 
Feature Extraction software, version A.7.5.1 (Agilent 
Technologies). The microarray data has been submitted 
to the GEO database (accession number GSE31605). 

Data analysis 

The microarray data were first subjected to linear nor- 
malization to allow comparison between arrays. All 
probes were compared to the reference sequence using 
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Figure 1 Experimental workflow and numbers of differentially expressed probes (DEPs) resulting from BMP4 or BMP7 treatment. (A) 

Seven breast cancer cell lines were cultured on 24-well plates, allowed to adhere for 24 h, and treated with the BMP ligand or vehicle for 30 
min, 1 h, 3 h, 6 h, 12 h, and 24 h. Experiments were performed in triplicate, and collected cells were pooled. (B) The expression data from each 
cell line were individually filtered according to the following criteria: differential expression of at least 2-fold at a minimum of one time point. 
The number of DEPs per cell line is represented. (C) The expression data from individual time points of every cell line were filtered according to 
a 2-fold cutoff in expression change. The number of DEPs per time point is shown. 



BLAST (v.2.2.23). Ensembl IDs for the probes were 
obtained by examining the probe's genomic location. 
The Ensembl Homo sapiens database version 60.37e was 
used. This process resulted in the annotation of 84% of 
all the probes. A total of 66% of the probes mapped 
uniquely to genes and 16% mapped to multiple genes. 

In order to determine differentially-expressed genes, 
expression data were subjected to three types of filter- 
ing: cell line-specific, time point-specific and general fil- 
tering. Cell line-specific filtering was done separately for 
each cell line, following the criteria of a differential 
expression of at least 2-fold in a minimum of one time 
point. In the time point-specific filtering, data from 
each time point were independently filtered according 
to a 2-fold expression change cutoff. General filtering 
was performed on all the data from all the cell lines 
together (separately for BMP4 and BMP7) following the 
next criteria: probes with a differential expression of at 
least 3-fold in at least three events and/or 2-fold in at 
least four events were considered for subsequent 



analysis. An event refers to any time point of any cell 
line, resulting in a maximum number of 30 events (5 
cell lines and 6 time points per cell line). The data sets 
produced by general filtering were further hand-anno- 
tated to reduce the number of probes with multiple 
annotations. Uniquely annotated probes are designated 
hereafter as "genes, " whereas the terms "probe" and 
"genetic element" refer to multiple annotated probes or 
any data including them. The gene lists resulting from 
general filtering were ranked according to the number 
of events in which they showed regulation. Furthermore, 
all the probes derived from general filtering were sub- 
jected to hierarchical clustering using correlation 
metrics, agglomerative strategy and average linkage 
method. 

Enrichments of gene ontology (GO) terms were per- 
formed on several data sets applying Fishers exact test 
and using all genes present on the microarray as a refer- 
ence [22]. In all the GO enrichment analyses, only 
probes with unique annotation were used. 
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A model-based clustering method [23] that allows the 
finding of clusters of genes with similar expression pro- 
files was performed using MCLUST R package [24]. 
This method was applied to data resulting from cell 
line-specific filtering. Data were log2-transformed and 
scaled to unit length. In the model-based clustering 
method, the clusters are considered to be groups dis- 
playing multivariate distributions. Several models were 
fitted to the data. The selections of the best model and 
the number of clusters were made based upon maximiz- 
ing Bayesian Information Criterion (BIC) values for the 
specific model and number of clusters that best repre- 
sented the data. The data analyses were performed 
using the Anduril data analysis framework [25] and R 
[26]. 

Results 

The aim of this study was to uncover the transcriptional 
responses of BMP4 and BMP7 signaling in breast can- 
cer. To this end, we selected breast cancer cell lines 
with low endogenous expression of BMP4 (HCC1419, 
SK-BR-3), BMP7 (MDA-MB-231, T-47D) or both 
(HCC1954, MDA-MB-361, ZR-75-30) [12,14] and trea- 
ted them with the corresponding BMP ligand (rhBMP4 
or rhBMP7) and vehicle controls (Figure 1A). Global 
gene expression levels were analyzed at six different 
time points from 30 min to 24 h in order to reveal the 
temporal patterns of transcriptional changes. 

Overall transcriptional response to BMP4 and BMP7 
treatment 

Due to the multidimensional nature of our data, we 
used three different filtering approaches, each of them 
allowing analysis from a different perspective. Cell line- 
specific analysis of the expression data evidenced con- 
siderable variation in the number of differentially 
expressed probes (DEPs) from one cell line to another 
(Figure IB), implicating distinct differences in their tran- 
scriptional response to BMPs. Further evaluation of 
these results revealed that BMP4 treatment resulted in 
greater amounts of DEPs than BMP7 (average number 
of DEPs per cell line: 5, 469 versus 3, 898 for BMP4 and 
BMP7, respectively; Figure IB). This finding could not 
be explained by the differences in cell lines used to 
study the two ligands, as a similar outcome was 
observed in the three lines treated with both BMP4 and 
BMP7 (HCC-1954, MDA-MB-361 and ZR-75-30). Time 
point-specific filtering revealed clear temporal variation 
in the number of DEPs (Figure 1 C). Generally, there 
was a tendency towards a greater amount of DEPs at 
later time points. In order to focus our attention espe- 
cially on those genes whose expression was most consis- 
tently and extensively affected by BMP4 and BMP7 
signaling, we performed a general filtering of the 



expression data according to the following criteria: fold 
change (FC) > +/-3 in at least 3 events and/or FC > 
+/-2 in at least 4 events. This resulted in the identifica- 
tion of 2, 421 and 1, 263 differentially expressed gene 
elements (1, 678 and 905 uniquely annotated probes) 
for BMP4 and BMP7 experiments respectively, further 
evidencing a more prominent effect of BMP4 than 
BMP7 on gene transcription. 

Unsupervised hierarchical clustering on the data sets 
resulting from general filtering (BMP4 and BMP7 sepa- 
rately) revealed that the samples originating from a par- 
ticular cell line mainly clustered together (Figure 2), 
suggesting considerable variation in the response of 
individual cell lines to BMPs. The most obvious exam- 
ples are MDA-MB-361, ZR-75-30, and HCC1419 (for 
BMP4) as well as MDA-MB-361 and HCC1954 (for 
BMP7). Clustering according to time point was an 
uncommon phenomenon, but it was observed for the 
samples derived from MDA-MB-231, ZR-75-30, 
HCC1954 and T-47D after 30 min of BMP7 treatment 
(Figure 2B). At the probe level, both BMP4 and BMP7 
hierarchical trees revealed a small subset of gene ele- 
ments that clustered tightly together (named clusters A 
and B hereafter, Figure 2, blue boxes). However, it is 
important to note that the expression of the genes in 
these clusters was not altered in a similar fashion in all 
the cell lines; rather, they were upregulated in some cell 
lines and downregulated in others, showing diverse tem- 
poral patterns. For both ligands, these clusters appeared 
to dictate the division of the samples into two major 
tree branches (Figure 2). 

In order to obtain a general view of the cellular func- 
tions regulated as a result of BMP4 and BMP7 signaling 
in breast cancer cell lines, we performed a GO enrich- 
ment analysis on the data sets resulting from general fil- 
tering. As might have been expected, functional 
categories related to regulation of transcription were 
among the most highly enriched for both ligands. Addi- 
tionally, genes involved in organ development were 
abundantly regulated as a result of stimulation with 
either ligand. In this regard, BMP4 seemed more often 
to regulate genes involved in skeletal system develop- 
ment (Table 1), while BMP7, on the other hand, 
appeared to regulate genes involved in epithelial devel- 
opment, neurogenesis and tube development (Table 2). 

A common synexpression group of genes regulated in 
response to BMP4 and BMP7 signaling 

As mentioned above, hierarchical clustering unveiled 
gene clusters A (BMP4, containing 329 probes) and B 
(BMP7, 228 probes) with highly correlated expression 
patterns (Figure 2, blue boxes). Interestingly, of all the 
probes contained in clusters A and B, 210 (154 genes 
with known and unique annotation, named group C 
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Figure 2 Unsupervised hierarchical clustering of microarray data. The data sets resulting from general filtering (2, 421 probes for (A) BMP4 
and 1, 263 probes for (B) BMP7) were subjected to this analysis. For both ligands, there is an evident cluster of gene elements showing very 
highly correlated expression patterns throughout the samples (gene clusters A and B, blue boxes). 



hereafter, Additional file 1) were present in both clusters 
and thus represent shared BMP target genes. Direct 
comparisons of the expression patterns of group C 
probes in the three cell lines (HCC1954, MDA-MB-361 



and ZR-75-30) treated with both ligands revealed ele- 
ments of similarity between BMP4 and BMP7 response 
in the same cell line but high variability between differ- 
ent cell lines (Figure 3). For example, these genes were 
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Table 1 Enriched gene ontology categories for differentially expressed genes as a result of BMP4 treatment. 

Category Number of genes p-Value 

Biological process 

GO:0006355: regulation of transcription, DNA-dependent 240 0.005 

GO:0009888: tissue development 128 0.014 

GO:0001501: skeletal system development 53 0.031 

GO:0030154: cell differentiation 237 0.041 
Molecular function 

GO:0043565: sequence-specific DNA binding 94 0.043 

GO:0003700: sequence-specific DNA binding transcription factor activity 133 0.046 



commonly upregulated in MDA-MB-361 while downre- 
gulated in HCC1954. The GO enrichment analysis 
unveiled 23 enriched biological process terms, of which 
21 could be classified into two functional categories, 
namely development and morphogenesis (16 terms) and 
gene expression (5 terms) (Additional file 2). 

Temporal patterns of transcriptional response to BMP 
signaling 

Model-based clustering analysis of the expression data 
was performed to distinguish clusters of genes with 
similar temporal profiles of expression. With this 
method we identified 12 to 22 probe clusters for each 
cell line and these clusters could be subsequently classi- 
fied into four main categories (Tables 3 and 4 and Addi- 
tional file 3). Gene elements that were first regulated at 
30 min or 1 h were classified as early, 3 h or 6 h early- 
intermediate, and 12 h late-intermediate responders, 
regardless of their expression at later time points. Late 
responders included those probes differentially 
expressed exclusively at the 24 h time point. Representa- 
tive examples of clusters in the different temporal cate- 
gories are depicted in Figure 4 A and 4B. 
We had already explored the biological functions of the 
genes differentially expressed upon BMP treatment. As 
we were now able to temporally classify the genes, we 



became interested in evaluating whether temporal pat- 
terns of expression and gene function could be related. 
Therefore, we grouped the clusters from each temporal 
category (early, early-intermediate, late-intermediate, 
and late) and performed GO enrichment analyses on the 
genes within these groups (Additional file 3). Finally, 
results from the different cell lines were combined. 
Although enriched GO terms were not found in every 
temporal stage of every cell line, many GO terms were 
enriched in the four temporal categories for BMP4 (Fig- 
ure 4 C). Of all the functional terms enriched through- 
out the experiments, those related to development were 
especially abundant in the early and early-intermediate 
phases. Terms connected with regulation of gene 
expression appeared at all times, although they were less 
abundant at the late stage of 24 h. Metabolism-asso- 
ciated terms were also present throughout the experi- 
ment, but were most prominent at late-intermediate 
time points. Signal transduction appeared to be an 
affected biological process in all but the late-intermedi- 
ate phases; it seemed most profoundly altered during 
early stages. Many terms related to cell proliferation and 
DNA repair were enriched among late responder genes, 
but it should be noted that all these terms emerged 
from one single cell line, MDA-MB-361. The results 
from BMP7 data were very limited and may be 



Table 2 Enriched gene ontology categories for differentially expressed genes as a result of BMP7 treatment. 



Category 



Number of genes 



p-Value 



Biological process 

GO:0030182: neuron differentiation 

GO:0048730: epidermis morphogenesis 

GO:0051239: regulation of multicellular organismal process 

GO:0045944: positive regulation of transcription from RNA polymerase II promoter 

GO:0030855: epithelial cell differentiation 

GO:0016481: negative regulation of transcription 

GO:0035295: tube development 

Molecular function 

GO:0043565: sequence-specific DNA binding 

GO:0003705: RNA polymerase II transcription factor activity, enhancer binding 
GO:0016564: transcription repressor activity 
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Figure 3 Hierarchical clustering heat map of group C genes with supervised clustering at the sample level. Data from each cell line are 
grouped, and time points are arranged in temporal order from left to right: 30 min, 1 h, 3 h, 6 h, 12 h and 24 h. 



summarized as an enrichment of genes involved in regu- 
lation of metabolism and gene expression among the 
early-intermediate and late time categories. 

Potential novel BMP target genes in breast cancer 

One of the main goals of this study was to identify new 
BMP4 and BMP7 target genes relevant in breast cancer. 
In order to distinguish those genes most often and 



ubiquitously regulated by BMP treatments, differentially 
expressed genes resulting from general filtering (1, 678 
and 905 for BMP4 and BMP7 experiments, respectively) 
were ranked according to the number of times a gene 
was up- and downregulated throughout the series of cell 
lines and time points (Tables 5 and 6 and Additional 
file 4). It is interesting that although the proportion of 
up- and down-regulation events was roughly equal 
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Table 3 Summary of the temporal clusters obtained from the analysis of BMP4 data. 
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when considering all the differentially expressed genes, 
induction of gene transcription clearly prevailed over 
inhibition when considering the 100 top-ranked genes 
(75% and 73% of the events for BMP4 and BMP7, 
respectively). For BMP4, 80 genes were regulated in 10 
or more events, while for BMP7, the analogous number 
of genes was 29 (Tables 5 and 6). Out of the 30 possible 
regulation events for a single gene (5 cell lines and 6 
time points), the actual maximums were 23 events 
(PTPRG; protein tyrosine phosphatase, receptor type, G) 
for BMP4 and 19 events (GNRHR; gonadotropin-releas- 
ing hormone receptor) for BMP7 (Tables 5 and 6); all of 
these were upregulation events. As expected, several 
members of the Id family of inhibitors of DNA binding, 
well-known targets of BMPs [27], were strongly induced 
by both BMP ligands. In addition to PTPRG and 
GNRHR, other genes strongly upregulated by both 



ligands included APOC2 (apolipoprotein C-II) and an 
as-yet unnamed gene encoding an uncharacterized pro- 
tein {C12orf42). DUSP2 (dual specificity phosphatase 2) 
and MAP3K5 (mitogen-activated protein kinase kinase 
kinase 5) were highly induced only by BMP4. BMP7, on 
the other hand, strongly promoted the expression of 
genes including PBX1 (pre-B-cell leukemia homeobox 1) 
and ZSCAN4 (zinc finger and SCAN domain containing 
4), which were not among the genes most intensely 
regulated by BMP4. 

Discussion 

In recent years, it has become increasingly accepted that 
the deregulation of mechanisms normally involved in 
developmental processes has tumorigenic effects in adult 
tissues. One example is the BMP family of growth factors, 
whose function in cancer physiology has been 
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Table 4 Summary of the temporal clusters obtained from the analysis of BMP7 data. 
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Gene clusters were classified in four temporal categories (early, early-intermediate, late-intermediate and late). Clusters with unclear profiles were classified as 


undetermined. The cluster ID column contains the ID number of the cluster. 







demonstrated in many tumor types, including breast can- 
cer [9-11]. In spite of this, little is known about BMP tar- 
get genes in the context of tumors. The transcriptional 
responses of breast cancer cells to BMP signaling have 
been studied only minimally. More precisely, the effects of 
BMP2 and BMP7 treatments on transcription in MCF-7 
and MDA-MB-468 breast cancer cell lines, respectively, 
have been analyzed using cDNA microarrays of limited 
content (from several hundreds to 14, 500 gene probes) 
[28-30]. In this study, we have therefore set up an experi- 
mental procedure to identify potential BMP target genes 
in breast cancer by studying the effects of two BMP 



ligands, BMP4 and BMP7, on genome-wide gene expres- 
sion. These two BMPs were selected based on their essen- 
tial role in breast cancer, which we and others have 
demonstrated in recent years [9-11]. Both ligands are 
highly expressed in primary breast carcinomas as well as 
in breast cancer cell lines [12-14,16,17]. BMP7 expression 
was also shown to be associated with early bone metastasis 
[15]. Additionally, in vitro studies have implicated BMP4 
and BMP7 as important regulators of proliferation and 
migration of breast cancer cells [14,18,20,21,31]. 

Our experimental approach allowed multiple types of 
analyses and revealed interesting insights into how the 
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Figure 4 Time series analyses. Representative clusters from the four temporal categories are shown for BMP4 (A) and BMP7 (B). For each 
cluster, the upper figure shows the levels of differential expression through the time series for every probe. The lower chart represents the 
average value of differential expression for all the probes in the cluster along the time scale. The number of probes in each cluster is indicated 
under the cell line name. (C) GO analyses of the four temporal categories were performed, and data from the five cell lines were combined. The 
enriched GO terms for BMP4 are depicted. 
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Table 5 Top most regulated genes after BMP4 signaling. 
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Table 5 Top most regulated genes after BMP4 signaling. (Continued) 
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ATP/GTP binding protein-like 1 
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G protein-coupled receptor 20 
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Differentially expressed genes resulting from the general filtering of BMP4 data (1, 678) were ranked according to the number of times a gene was up- and 
downregulated throughout the series of cell lines and time points. Only those genes with total rank value of at least 10 are listed. 



stimulation of BMP4 and BMP7 signaling influences the 
transcriptome of breast cancer cell lines. First of all, we 
showed relatively high numbers of DEPs as a result of 
BMP treatments, indicating a strong impact of BMP4 
and BMP7 on the cell lines studied. Moreover, the tran- 
scriptional response to BMP4 was of a clearly higher 
magnitude than that induced by BMP7. Another aspect 
of the study was the opportunity to compare the effects 
of BMP signaling between cell lines. Interestingly, clear 
differences were seen in the amounts of DEPs, as well as 
in their expression patterns, as revealed by hierarchical 
clustering. BMP signaling pathways are regulated in a 
very complex manner and at many different levels, from 
the availability of BMP receptors, BMP ligands and 
BMP antagonists in the extracellular compartment to 
the presence or absence of various intracellular signal 
mediators and transcriptional co-activators or co-repres- 
sors [32,33]. Therefore, multiple factors influence the 
outcome of BMP signaling on the transcriptional level 



in a given cell. We have previously reported that all six 
BMP specific receptors (ACVR1, BMPR1A, BMPR1B, 
ACVR2A, ACVR2B, and BMPR2) are uniformly 
expressed among the breast cancer cell lines studied 
here [12]. Similarly, we have shown that SMAD4 is 
expressed and that phosphorylation of SMAD- 1/5/8 is 
induced in these cell lines after BMP7 and BMP4 treat- 
ment [14,21], Taken together, the expression profiles of 
BMP specific receptors or the mediators of the canoni- 
cal intracellular pathway do not seem to have a major 
role in explaining the different transcriptional responses 
in the breast cancer cells. Nevertheless, due to the com- 
plexity of BMP signaling regulation, it is easy to under- 
stand that different cell lines may have different 
transcriptional responses to BMP stimulation. This 
observation highlights the importance of testing multiple 
cell lines when studying BMP signaling in cancer. An 
additional finding that could be inferred from our data 
is that induction of gene transcription, compared with 
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Table 6 Top most regulated genes after BMP7 signaling. 

# regulated events 
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10 


ENSG00000164104 


HMGB2 


high-mobility group box 2 


8 


2 


10 


ENSG00000250589 


DUX4 


double homeobox 4 


2 


8 


10 


ENSG000001 30559 


CAMSAP1 


calmodulin regulated spectrin-associated protein 1 


3 


7 


10 


ENSG000001 66501 


PRKCB 


protein kinase C, beta 


4 


6 


10 


ENSG00000186115 


CYP4F2 


cytochrome P450, family 4, subfamily F, polypeptide 2 


4 


6 


10 



Differentially expressed genes resulting from the general filtering of BMP7 data (905) were ranked according to the number of times a gene was up- and 
downregulated throughout the series of cell lines and time points. Only those genes with total rank value of at least 10 are listed. 



inhibition, was the common response among those 
genes most frequently regulated by BMP4 and BMP7 in 
breast cancer cells. Likewise, previous microarray-based 
transcriptomic analyses of TGF-P and BMP have shown 
that induction of gene expression is the predominant 
response of mammalian cells to stimulation by these 
growth factors [28,30,33-35]. 

After BMP4 and BMP7 stimulation, the microarray 
analyses identified a large number of differentially 
expressed genes in our panel of cell lines. To explore 
the biological functions of these genes, GO enrichment 
analyses were performed. These revealed very similar 
results for both BMP ligands, namely, regulation of tran- 
scription and developmental processes. It seems, there- 
fore, that the functions most prevalently influenced by 
BMP signaling in breast cancer cells do not differ 
remarkably from conventional roles that BMPs possess 
during development [36,37]. 



Synexpression groups are synchronously coexpressed 
gene sets, particularly apparent during embryonic devel- 
opment and in the response of cells to hormones and 
growth factors [38,39]. Our analyses unveiled that treat- 
ment of breast cancer cells with either BMP4 or BMP7 
resulted in the coordinated expression of a group of 
genes (clusters A and B, respectively). Most interest- 
ingly, a considerable number of the genes in these two 
synexpression groups were common for the two ligands 
(group C). Moreover, our data indicated that treatment 
of a cell line with either BMP4 or BMP7 results in simi- 
lar transcriptional responses of group C genes. We 
therefore hypothesize that group C represents molecular 
responses shared by the BMP4 and BMP7 signaling 
pathways. This finding prompted us to ask what func- 
tions these common genes fulfill in the cell. GO enrich- 
ment analysis of the genes in group C revealed that 
these genes are involved in two main biological 
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processes, regulation of gene expression and regulation 
of development and morphogenesis. These results sup- 
port the notion that genes known to regulate develop- 
ment also have functions that are important for the 
maintenance of cancer cells. 

We also studied the temporal patterns of the tran- 
scriptional response after BMP treatment. The number 
of DEPs showed a tendency to increase with time, a 
trend previously noticed in transcriptome analysis of 
TGF-P family members in murine mammary epithelial 
cells and in breast cancer cells [29,35]. The DEPs could 
be grouped according to their temporal pattern of 
expression, varying from early to late responders. These 
temporal clusters were found in every cell line, and 
some of them even contained over a thousand gene ele- 
ments. The next logical step was to explore whether 
there was a time-dependent shift in the distribution of 
gene functions. Although GO enrichment results were 
not obtained for all the probe clusters of all the cell 
lines, interesting features could be identified, especially 
in the case of the BMP4 data. Transcriptional regulation 
in the first 6 hours concentrated most notably on genes 
involved in developmental processes, metabolic pro- 
cesses, gene expression and signal transduction. Gene 
expression was also well-represented after 12 hours, 
while metabolism became by far the most prominent 
function at this time point. Most interestingly, 24 hours 
after BMP4 stimulation there was an evident overrepre- 
sentation of genes involved in cell proliferation, 
although this phenomenon was observed exclusively in 
MDA-MB-361 cells. All in all, the enriched biological 
functions indeed fluctuated in time and in a logical 
sequence, with regulation of gene expression and signal 
transduction leading to changes in metabolism and 
finally to regulation of cell proliferation, a phenotype 
relevant for cancer cell physiology. The fact that we did 
not see enrichment of cell proliferation-associated func- 
tions in more than one cell line could be due to differ- 
ences in the speed of BMP signaling in different cell 
lines. Even though a longer experiment certainly could 
have clarified this issue, we concentrated our analysis on 
the first 24 hours after BMP treatment because we were 
interested primarily in the identification of BMP target 
genes. 

As mentioned, one of the main goals of this study was 
to identify potential novel gene targets of BMP signaling 
relevant in breast cancer. We provided lists of candidate 
genes that are strongly and rather uniformly regulated 
by BMP4 or BMP7 throughout the cell lines and time 
points. Some of them, such as members of the Id family 
of inhibitors of DNA binding, are well-known BMP tar- 
get genes [27,29]. Id proteins are transcription factors 
that regulate cell growth and differentiation [40], and all 
four members of the protein family play crucial roles in 



various aspects of normal and malignant breast biology 
[41]. Others are newly linked, in this work, to BMP sig- 
naling, and some of these genes have interesting con- 
nections with tumor biology, such as PTPRG or DUSP2. 
A positive feedback regulation where BMP treatment 
leads to increased expression of BMP antagonists is 
known to exist. In our study, no consistent expression 
changes were observed for any of the known BMP 
antagonists, such as noggin, gremlin, sclerostin and fol- 
listatin. Previous studies have shown a wide time win- 
dow in the induction of e.g. noggin expression in 
different tissues after BMP treatment, ranging from 1 to 
48 hours [42,43]. Thus it is possible that the feedback 
effect in the breast cancer cells was not evident at time 
points analyzed here. 

Protein tyrosine phosphatases (PTPs) are key regula- 
tors of the cellular protein phosphorylation balance, cri- 
tical in the control of a wide spectrum of physiological 
processes such as cell proliferation, differentiation, 
transformation, transport and locomotion. Subsequently, 
aberrations in phosphorylation processes play a major 
role in the pathogenesis of numerous diseases, including 
cancer [44,45]. PTPRG is a receptor-type PTP impli- 
cated as a candidate tumor suppressor gene in several 
types of tumors, including breast cancer [46,47]. In 
MCF-7 breast cancer cells, PTPRG inhibits proliferation 
and anchorage-independent growth and reduces tumor 
formation in a xenograft model [47,48]. Delayed cell 
cycle re-entry by increasing the level of cell cycle regula- 
tors p21 and p27 could explain the inhibitory effect of 
PTPRG on cell growth [47]. Based on the above, upre- 
gulation of PTPRG in BMP-stimulated cancer cells 
could contribute to the observed BMP-induced antipro- 
liferative effect [14,21]. DUSP2 also belongs to the PTP 
family of phosphatases. It is a mitogen-activated protein 
kinase (MAPK) phosphatase (MKP) that dephosphory- 
lates both threonine and tyrosine residues within target 
MAPKs leading to their deactivation. MAPK signaling 
controls cellular processes such as proliferation, differ- 
entiation, migration and apoptosis [48]. Therefore, 
abnormal MKP activity, and hence anomalous MAPK 
signaling, has important consequences for processes cri- 
tical to the development and progression of human can- 
cer. The role of DUSP2 in cancer has been examined in 
only a few studies, and data are controversial. Overex- 
pression of DUSP2 expression was found in 37 of 39 
malignant effusions from serous ovarian carcinoma 
patients and was associated with poor survival [49]. By 
contrast, decreased DUSP2 transcript levels were 
reported in cancerous breast, colon, lung, ovary, kidney 
and prostate tissues, and reduced DUSP2 protein levels 
were observed in cervical and colon cancer [50]. Addi- 
tionally, DUSP2 suppression was associated with tumori- 
genesis and malignancy in colon cancer, and DUSP2 
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overexpression induced apoptosis and inhibited tumor 
growth in HeLa cells in vitro and in xenograft models 
[51]. 

Conclusions 

All in all, we show that BMP4 and BMP7 have strong 
effects on gene expression in breast cancer cells, and 
that this transcriptional response and its functional out- 
come follow a temporal sequence. Our data support the 
existence of a synexpression group of regulated genes 
that represent the core molecular responses shared by 
BMP4 and BMP7 signaling pathways. Additionally, we 
provide a list of potential novel BMP target genes rele- 
vant in breast cancer. 

Additional material 



Additional file 1: List of probes in group C. The expression values are 
represented as log2 ratios. Data was colored according to the 
accompanying color key. The data sheet contains 210 probes. 

Additional file 2: GO enrichment results for group C genes. 

Additional file 3: Model-based clustering results and GO 
enrichment. The probe clusters obtained are represented as well as the 
results of the GO enrichment analyses. GO terms are colored according 
to the accompanying color key. 

Additional file 4: Most commonly regulated genes with rank 6 or 
higher. The expression values are represented as log2 ratios. Data was 
colored according to the accompanying color key. The data sheets 
contain 608 and 281 genes for BMP4 and BMP7, respectively. 
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